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Abstract 

We construct a new class of non-topological solitons in renormalizable scalar 
field theories with non-linear self-interactions. For large charge Q, the soli- 
ton mass increases linearly with Q, t.e., the soliton mass density is approxi- 
mately independent of charge. Such objects could be naturally produced in a 
phase transition in the early universe or in the decay of superconducting cosmic 
strings. 
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Non-topological solitons (hereafter, NTS’s) are solutions of classical field theories which 
are stable by virtue of a conserved Noether charge carried by fields confined to a finite 
region of space 1-6 . Recently, these solutions have been studied under the guise of Q 
balls 7 , cosmic neutrino balls 8 , quark nuggets 9 , and soliton stars 10 , and a scenario for 
producing them in a phase transition in the early universe has been suggested 11 . The 
simplest renormalizable theory with NTS solutions is an unbroken global U(l) theory of 
two coupled scalar fields. In this paper, we study the classical NTS solutions of this theory 
in detail. We will confine our discussion to NTS masses M C rn^i/m 2, (where m is the 
mass of the scalar fields) , so that gravity can be neglected. The inclusion of gravitational 
effects for large mass solutions has been studied elsewhere by one of us 12 . 

Consider the Lagrangian for a real scalr field a and a complex scalar $, 

£ = i(a^) 2 + |^*| 2 -C/(|*|,<7). (1) 

The most general interaction invariant under the discrete symmetry a — * —a and the 
global 1/(1) symmetry $ — ► e*“$ is 

£'(1*1.0) = j(o 2 - + 9 2 |$| 4 + ?%|*| V (2) 

* °o 

We are interested in the case Oq > 0, so that the discrete symmetry is spontaneously 
broken in the ground state, < a >= ±<tq. Although this leads to the familiar domain 
wall problem 13 (if oq > 1 MeV), we assume this can be circumvented, e.g., by adding 
explicit discrete symmetry-breaking terms 14 . (Since these terms can be very small, we 
neglect them in our analysis.) Alternatively, we may consider a as the real component of 
a complex field whose vacuum expectation value breaks a global C/(l) / symmetry; such a 
theory has no stable domain walls, but may have stable vortices (cosmic strings). In either 
case, the 1/(1) symmetry carried by $ is unbroken, so there is a conserved current 
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^ and an associated conserved charge 

Q = J StJq. (4)* 

(The conserved charge stabilizes the NTS against decay.) In vacuum, the particle masses 
are given by 

ml = 2A<7q, + m\ 0 . (5) 

Spherically symmetric NTS solutions for the theory of Eqns. (1-2) were first studied 
by Friedberg, Lee and Sirlin 3 for the special case = g = 0. For large Q, they are 

characterized by an interior “false vacuum” region where a cz 0, surrounded by a thin 
domain wall where a rapidly approaches its ground state value a = oq. In the NTS 
interior, the potential energy density in a is balanced by the pressure of the massless $ 
charges, which are confined by the mass gap ~ at the domain wall. 

Here, we consider a more general class of solutions. In particular, we show below 
that inclusion of the repulsive j$| 4 self- interaction completely changes the behavior of 
the NTS’s, altering the mass-radius and mass-charge relations. In the next Section we 
discuss approximate analytic NTS solutions, valid for large Q. In Section III, we construct 
numerical solutions, which are useful for displaying the behavior at small charge. We 
consider the spectrum of NTS excitations (phonons and particles) in section IV. In the 
Conclusion, we speculate on the place of NTS’s in particle physics models and briefly 
outline how such objects might be produced in the early universe. 

II. Analytic NTS Solutions 

From Eqns. (3-4), a classical solution $(r, t) must be time-dependent for the system 
to have non-zero charge. The lowest energy solutions are of the spherically symmetric 
form 1-4 

= a(r,t)=o(r). (6) 
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It proves convenient to define the rescaled field variables, 




spatial coordinate, 


y = ga 0 f ; 


( 7 ) 

( 8 ) 


masses 


fho _ fh$ a tZ> _ 1 

m a Tn $0 w goQ ’ 


and coupling constant, 



The energy functional for the solution (6) can then be written 


( 9 ) 


( 10 ) 


M = 


Mao 47TCQ 

9 ~ 9 


y 2 dyfyd y <p) 2 + ^{d y a) 2 + ^w 2 <p 2 + U[fi > o)} > ( 11 ) 


where 


U{<P,o) = ^(m|* + m\ a b 2 )v 2 + ^- + ^(a 2 - l) 2 


(lla) 


is the rescaled potential. The soliton charge is 



( 12 ) 


The rescaled field equations are 

^ + y^y + = ^C**"* 2 + ^ 2 ) ( 13 ) 

d 2 d 2 do r _._o , .0 .2- / X 

V + ^ - A<, <" - = m *°' e < H > 


The method of constructing NTS solutions is, in principle, straightforward: they are 
uniquely determined by specifying the parameters w, m$$, and A, and by imposing 
the boundary conditions v?(y) — ♦ 0, d(y) — * 1 as y — *• oo, and d<p/dy,da/dy — ♦ 0 as y — ► 0. 
One then solves Eqns. (13-14) subject to the constraint that M, the rescaled NTS mass, is 


4 



ORIGINAL page is 
OF POOR QUALITY 


minimized; this also fixes the soliton charge Q. We note that each such solution actually 
yields a two-parameter family of NTS’s, since g and <to have been scaled out. 

To gain insight into the nature of the NTS solutions, one can think of Eqns. (13) 
and (14) as the equations of motion of a particle rolling in a two-dimentional potential 
V e ff(<p,a), with “space coordinates” (<p, d) and “time coordinate” y. The particle is 
subject to a friction force which varies as y -1 . The effective potential for this mechanical 
analogy is 

= ~U{<p,o) (15) 

and is shown in Fig. 1 for a particular choice of parameters. At y = 0, the NTS center, the 
particle starts at rest, while at late “times”, y — > oo (».e., far outside the NTS), it must 
approach the vacuum state d = l,<p — 0 on the ridge of the potential. Since the particle 
rolls with friction, it must start out at a point higher on the potential saddle, where d < 1, 
<P> o. Clearly, then, a necessary condition for the existence of a solution is 

maxV e ff(<p,a) > V e f f { 0,1) = 0. (16) 

Suppose this condition is true for some choice of the parameters, and then imagine 
decreasing w while keeping the masses and couplings fixed. From Eqn. (15), the inequality 
(16) then approaches equality so that, for a solution to exist, the friction term must 
become negligible. But approximate solutions in which the friction term can be consistently 
dropped are easy to construct. They are of the form 

«*>“{£ y y >Y <«.) 

*M“{? / y>Y. W 

Here, the particle stays close to its initial position (<p c ,d c ) until a large time, y — Y, when 
it suddenly rolls down the potential and asymptotically comes to rest at the vacuum state. 
For the ansatz (17), the friction terms y~ l d<p/dy, y~ 1 do/dy would only be appreciable 
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where the field gradients are large, at y ~ Y; however, if Y is sufficiently large, these 
terms are suppressed. 

The solution above describes large, thin-walled solitons of physical radius R = Y / goa, 
with negligible surface energy. We call these large Q solitons “bag” solutions by analogy 
with hadron models. One feature of these bag solutions should be emphasized: by con- 
struction [Eqns. (8-10)], they only exist for g ± 0, *.e., for non-zero |$| 4 coupling. To see 
this, suppose we instead define a different coordinate rescaling x = oqt, and correspond- 
ingly define rescaled parameters m$$, w , and A, obtained by setting g = 1 in Eqns. 

(9) and (10). Then Eqn. (13) becomes 

“f + + (™ 2 ~ = <p{m\ a b 2 + g 2 & 2 ). (18) 

In the limit g — ♦ 0, the solution (17a) becomes inconsistent. Rather, for g = 0 and 
a = b c = const. Eqn. (18) is just the linear equation for spherical waves, with solution 3 

<P = ~ sin [(to 2 - - mj^dj) 1 / 2 *] . (19) 

For large Q, it thus turns out that, unless g is very small, g C Q -1 / 4 , the solutions with 
non-zero g are qualitatively different from those with g = 0. In particular, at large Q, for 
g = 0, the central amplitude <p c grows as a power of Q while, for g 0, <p c approaches 
a constant, independent of Q (see below). This arises because the <7 2 |$| 4 contribution to 
the NTS energy density suppresses the amplitude of $. We will compare the cases <7 = 0, 
g 7^ 0 more fully in our numerical work below. 

Inserting the ansatz (17) in Eqns. (11) and (12), we find the rescaled charge 

Q = yV’uite?, (20) 

and NTS mass 

M = Yy 3 {^tZ 2 £ 2 + U{<p e ,be)}, (21) 
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Using Eqn. (20) to eliminate w, Eqn. (21) can be rewritten 




( 22 ) 


where V = ^pF 3 is the rescaled NTS volume. Minimizing the energy with respect to V at 
fixed charge gives 


V = 


M = 


Q\J 2U(£ c ,d c ) 

Vc 


<p\j2U{<p c ,Oc) 

If we minimize the mass M (at fixed Q) with respect to <p c and cr c » we obtain 


(23) 


&c =0, ip c = A 1 / 4 , 


(24) 


so that <p\ = Aoq jg 2 . Thus, the NTS energy is 


M=g(A 1 / 2 + m | $ ) 1 /2 = t5o Q 

(25) 

or 


M = w 0 Q = Q{\ l / 2 gal + m^) 1 / 2 . 

(26) 

For Eqn. (23), the NTS charge is 


Q - w 0 <p 2 V, 

(27) 

so the mass-radius relation is 


\l/2-2 m 2 

9 

(28) 


i.e., large charge NTS’s have approximately constant density, p = Wq<p 2 . We note that 
these relations, i.e., M ~ Q,R ~ Q 1 / 3 , are very different from the “free” case g = = 0 

treated in Ref. 3, where M ~ Q 3 / 4 , R ~ Q 1 / 4 . In particular, the solutions displayed above 
have a minimum frequency, w — ► iwq as Q — ► oo, while “free” solitons have w = n/R ~ 
q-1/4 Yqx the self-interacting theory, the NTS solutions are very akin to Q-balls 7 . 
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Non-topological solitons are quantum-mechanically stable if they are the lowest energy 
configuration of fixed charge. Using Eqn. (5), the stability conditions are 


M{Q) < Q(m|$ + rnjj 1 / 2 


and 


d 2 M 

dQ 2 


< 0 . 


(29) 


(30) 


Eqn. (29) ensures stability against decay into free $ particles, while Eqn. (30) expresses 
stability against fission into smaller soliton fragments. 

For large Q, the stability condition (29) becomes (from (26)] 

g\ l/2 ol < (31) 


independent of the “bare” mass To check stability against fission, we would need 

to include the surface terms neglected above. This is easily done: treated perturbatively, 
the surface energy 4irR 2 p (where p is the surface energy per unit area) contributes a term 
~ Q 2 / 3 to the NTS mass, and therefore d 2 M/dQ 2 < 0. Thus, the binding energy per unit 
charge increases monotonically with Q, and large Q NTS’s do not fission spontaneously. 
The addition of the surface energy also implies that the stability bound (29) is violated 
at sufficiently small charge, there exists a minimum charge Q m in above which NTS’s 
are absolutely stable. Generally, for Q ~ Q m ia the surface terms cannot be treated as 
perturbation, so we defer discussion of this point to the numerical solutions below. 

III. Numerical NTS Solutions 

To construct solutions at arbitrary Q, we have numerically solved Eqns. (13) and (14) 
subject to the boundary conditions noted in the previous Section. The program uses the 
parameters tu, rh$ a and A as inputs and returns M and Q as outputs. It is useful 

to write the energy functional as 

M — Ef( + Ep + E} y, (32) 


8 



where the gradient energy 


ORIGINAL, PAGE IS 

OF POOR QUALITY 


E k = Ait J y : 


d V \\{ d yt) 2 + \{ d y°) 2 


the “potential” energy 


Ep 


= Air J y 2 dy 


u(<e,5)~ 


w 2 <p 2 ‘ 


and the “number” energy 


Eft = 4?r J y 2 dyw 2 <p 2 = wQ. 


(Z3f 


(34) 


(35) 


For the analytic solutions of the previous Section, M a Efj. More generally, one can 
show that the NTS solutions must satisfy a virial theorem, 3 



(36) 


This relation is used to check the accuracy of the numerical solutions: the program iterates 
from a trial solution until the quantity 


Ep + \E_k 

Ep 


is sufficiently small, typically less than 10 4 or so. 


(37) 


In Fig. 2, we show two NTS solutions with parameters — 0, = 1, A = 

1/9 (assuming g ^ 0). The solid curves show <p(y), o(y) for w 2 = 0.38, with charge 
Q = 3.1 x 10 4 and mass M = 2.0 x 10 4 ; in this case, Wq = A 1 / 2 = 0.33, so this is a bag 
solution with w ~ u)q and M ~ wQ. The dashed curves show the solution for w 2 = 0.8, 
which has a charge Q = 74.5 and mass M = 72.2. For this case, the potential and surface 
energies are important and, as one sees from the figure, the solution is not at all bag-like. 
In particular, here d(0) = 0.5 instead of zero. As w is further increased, eventually we 
reach a charge Q m iu where the stability condition (29) breaks down. For the parameters 
above, this occurs for w 2 = 0.92, with <5 m in = 45.1 and M m ^ = 45.0. 
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It is interesting to consider what happens for Q < Q m ;„ . Although such NTS’s have 
positive binding energy, and are thus quantum mechanically unstable, classical solutions 
still exist. As xv 2 is increased from 0.92 (its value at Q m in), the charge and mass first 
continue to drop, but then M(Q) reaches a cusp at Q cr j t = 39.4, M cr j t = 39.5, u> 2 rit = 0.97. 
For w > ui cr it, the mass and charge grow with increasing w and, as ui approaches rh$ ff = 1 
from below, the NTS energy approaches the free particle energy Mf ree = m$ a Q = Q. 
Thus, a plot of M(Q) would show two branches, the upper one with ti; > u> cr i t , the lower 
with xv < u> cr it, which join at Q min . (Similar behavior at small Q was found in Ref 3.) 
Presumably, the upper branch solutions decay either to the lower branch or to free $ 
particles. 

For comparison, in Fig. 3 we show NTS solutions with the same parameters as above, 
but with <7 = 0, i.e., we have m$$ = 0, m$ a = m^/oo = 1, A = A = 1/9, and we plot 
<p(x), &(x) where z = oqt [see discussion around Eqn. (18)]. The solid curves are for the 
case u) 2 = w 2 /oq = 0.038, Q = Q — 9.6 x 10 3 , M = M/oq — 2.5 x 10 3 , while the dashed 
curves are for w 2 = 0.6, Q = 64.1, M = 59.4. It is clear that these solutions are far from 
bags, but are instead well approximated by the solutions (17b), (19) for large Q. For these 
parameters, the minimum charge for stability is = 34.5, M m j n = 34.5, which occurs 
at w 2 = 0.875. To see the effect of the <7 2 |$| 4 term graphically, we plot the soliton mass 
M(M) as a function of charge Q(Q) for the case g ± 0(g = 0), in Fig. 4. To compare the 
curves, we can take <7 = 1 for the model with non-zero g, so that we are plotting M/oq vs. 
Q in both cases. As expected, at large charge Q, the |$| 4 term increases the NTS mass. 
(We have displayed only the lower branch solutions.) 

To summarize, for the self-interacting theory, stable NTS solutions exist over the fre- 
quency range 

(A l / 2 gol + mis,) 1 / 2 < w < (m| CT + mis,) 1 / 2 (38) 

and for charges greater than a minimum charge Q m iu which depends on coupling constants. 
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IV. Soliton Excitations 

The phenomenology of non-topological solitons is determined in part by the small 
fluctuations about the NTS ground state. For example, at finite temperature, these fluc- 
tuations make an important contribution to the free energy of the NTS and are therefore 
crucial in determining the phase diagram of the theory. These excitations are of three 
types: i) a particle and $ particle-antiparticle excitations, with energies proportional to 
the (<r, $) masses in the NTS interior; ii) sound waves (phonons), with energies proportional 
to R~ iii) surface waves. In this Section, we compute the spectrum of NTS fluctuations 
of types (i) and (ii) for large Q. For simplicity, we shall focus on the case = 0. 

To study small amplitude fluctuations, we define the perturbed variables, 


-a 


a = a(r) + 6 a(f,t), 

1 / cos wt —sin ut 


sinwt 


cosut 


<p, 


(39) 

(40) 


where 


<p = 


(41) 


<p{r) + ^i(r,t)\ 

&<P 2 (r,t) J 

Here, o(r) and <p(r ) are the unperturbed NTS solutions of Eqn. (6). Linearizing in 
6a, 6<pi,6<P2 and using the unperturbed equations of motion for the background solution 
yields 

^ Jp “ + °2^ + *3b<Pi = 0 ( 42 ) 

0.360 + ^ Jp - V 2 + ai + <14^ 6<pi - 2w = 0 ( 43 ) 

d ( \ 

2u} dt 6,pi + ( dt* " = 0 . ( 44 ) 

Here, u> is the frequency of the background soliton, and 04, ...,<14 are functions of the 
unperturbed solution, 


„ _ 2 , 2 2 2 
a l — — 2~° + 9 <p “ w 


(«) 
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a 2 = *2 V 2 + M^ 2 °l) + 2A<7 2 

a o 

(46) 

2m la 

03 = 2 

° 0 

(47) 

04 = 2 g 2 <p 2 . 

(48) 


For arbitrary charge, one can in principle solve Eqns. (42-44) for the normal modes 
of the NTS; the solutions involve spherical Bessel functions and Legendre polynomials, 
with suitable boundary conditions imposed at the NTS surface. We shall instead consider 
modes of the large Q bags of Section II. 

Using the bag solution (24), we have 01 = 03 = 0, and 


TO^A 1 / 2 2 

a 2 — 

9 

(49) 

<14 = 2 u>q = 20A 1 // 2 <7q. 

(50) 


Thus, in the large Q limit, the So and 6<p modes decouple. Furthermore, since we are 
considering the limit of infinite NTS radius, we can approximate the normal modes with 
plane waves of the form e *(*° t ~ kx ). For the Sip modes, we then find a dispersion relation 
of the form 

fcg = A 2 + 3 wq ± u»o(4fc 2 + 9 wo) 1//2 > (51) 

where fc 2 = |£| 2 . At low momentum, i.e., k C u>o> the lower (acoustic) branch gives 
k,Q = it 2 / 3, which corresponds to the sound speed of a relativistic medium, c 2 = 1/3. The 
upper (optical) branch corresponds to $ particle-antiparticle excitations, with rest energy 
ko — y/&u>Q. 

For the So modes, in the bag limit the mode equation (42) is just the free Klein- 
Gordon equation with mass m = yfai [see Eqn. (49)]. We note that the stability condition 
(31) implies a 2 > 0 ; in other words, as one would expect, NTS stability corresponds to 
non-tachyonic So fluctuations. 
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Although we have studied the simplest theory containing non- topological solitons, 
many features of the solutions discussed here are universal. For example, for large charge, 
we found that scalar NTS’s have approximately constant mass density, i.e., mass propor- 
tional to their charge. This scaling holds whenever the charged field, be it fermion or 
boson, has a reduced mass in a region of non- zero vacuum energy. Thus, the linear scaling 
of mass with charge is valid for Q-balls, strange nuggets, bag models of hadrons, etc. 

Non-topological solitons can arise in any particle theory in which the lightest field 
carrying a global additively conserved quantum number gets a mass via a Higgs-type 
mechanism. Many extensions of the standard model, such as technicolor, grand unified 
theories, left-right symmetric theories and majoron models may have such a structure. 
Candidates for the requisite charged field include the lightest technibaryon and a massive 
neutrino. Whether the NTS solutions in such theories are actually stable depends only on 
ratios of coupling constants. Generally, one finds stability above some minimum charge 
Qmin (if & ii) > f° r theories with small coupling constants (of order unity) , <2 m in typically 
lies in the range between one and 10 4 or so. For Q m ^ < 1, we have the intriguing possibility 
that NTS’s jure the lowest lying solutions in the particle spectrum. 

In passing, let us note that NTS’s may also arise in theories with a multiplicatively 
conserved charge. An interesting example is supersymmetry, where Q is identified with R- 
parity. If the lightest supersymmetric particle is a scalar neutrino (or a Majorana fermion 
such as the photino), the associated NTS’s would be qualitatively similar to those described 
here. However, since R-parity is multiplicative, sneutrino solitons with R > 2 can decay 
by self-annihilation. Supersymmetric NTS’s would only survive if the NTS frequency u> is 
below the mass threshold of the lightest LSP annihilation channel or if i2mi n < 1. 

How might such objects come into existence? During a phase transition in the early 
universe, NTS’s may condense due to charge fluctuations or due to a charge asymmetry 
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(i.e., if the universe carries a net Q). 11 To study this process in detail, we are at present 
constructing NTS solutions at finite temperature. 15 Once formed, however, NTS’s are gen- 
erally vulnerable to evaporation 16 above a critical temperature at which they become the 
lowest free energy state. Whether they survive this evaporation phase is model-dependent. 

Another mechanism for NTS formation is the decay of cosmic strings. Suppose a is 
the real component of a complex field in a t/(l) x U(l)' theory. For certain choices of 
parameters, this theory is known to have vortex solutions which carry both current and 
charge. 17 In some cases, the resulting ‘vortons’ may evolve to a state in which the charged 
condensate is separated from the current; when this happens, it is likely that a stable NTS 
is formed. Since cosmic strings can be efficiently produced in phase transitions, 18 NTS’s 
may be produced in abundance. 
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Figure Captions 

1. Effective potential for the mechanical analogy, V e fj(<p,o) [Eqn.(15)]. 

2. NTS solutions <£(y),d(y), for = 0, rh$ a = 1, A = 1/9, j ^ 0. Solid curves: 

u 2 = 0.38, Q = 3.1 X 10 4 . Dashed curves: Cj 2 = 0.8, Q = 74.5. 

3. NTS solutions <p(x),a(x), for same parameters as Fig. 2, but with g = 0. Solid curves: 
u> 2 = 0.038, Q = 9.6 x 10 3 . Dashed curves: Q 2 = 0.6, Q = 64.1. 

4. NTS mass M/oq as a function of charge Q. Dashed curve: g = 0. Solid curve: g — 1 . 
Dotted curve: free particles, M = Qm$ a . 
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